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Ch ! Abstract 
i— » 

' We present the many-particle Hamiltonian model of Lipkin, Meshkov and 

Glick in the context of deformed polynomial algebras and show that its exact 
solutions can be easily and naturally obtained within this formalism. The 
(T) | Hamiltonian matrix of each j multiplet can be split into two submatrices 

associated to two distinct irreps of the deformed algebra. Their invariant 
\q • subspaces correspond to even and odd numbers of particle-hole excitations. 

o ' 
o 

^ . I. INTRODUCTION 

-i— > ■ 

In the sixties much interest has been devoted to formalisms for treating multiparticle 
systems and the quality of the approximations involved. To test the validity of the approxi- 



mations, quasi-exactly solvable models have been proposed (for a definition of a quasi-exactly 
solvable model see e.g. Ref. The comparison between the exact solutions and an ap- 
proximation could give a clear estimate of the quality of the approximation, which could 
further be applied to more complicated Hamiltonians. Among them of particular interest is 
the model of Lipkin, Meshkov and Glick (LMG) ||. Although simple enough to be solved 
exactly, in some cases the model is not trivial. Few analytic solutions (for a number of par- 
ticles up to 8) have been provided by LMG. Numerical solutions were also given for a larger 
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number of particles in the case where the total angular momentum reaches its maximum 
value. Here we study the exact solutions of the LMG model in the framework of deformed 
algebras. We first show that the LMG model corresponds to a deformed algebra of poly- 
nomial type and then we search for possible solutions associated with the representations 
of the corresponding deformed algebra. We show that the polynomial algebra introduces a 
new symmetry in the system, not known before, and this splits into two submatrices any 
Hamiltonian matrix to be diagonalized for a specific value of the angular momentum. More- 
over, it introduces a new quantum number which naturally distinguishes between an even 
and an odd number of particles. 

In the next section we recall the LMG model. In Sec. Ill we shortly introduce deformed 
algebras and describe the way one can use them to find exact solutions of the LMG model for 
an arbitrary given number of particles and for any specific value j of the angular momentum. 
In Sec. IV we presents algebraic and numerical solutions for a system with an odd (N = 
7) and an even (N = 8) number of particles. A general description of some supplementary 
solutions, inherent to the deformed algebra, is given in Sec. V. The last section is devoted 
to a summary. 

II. THE LMG MODEL 

As mentioned in the introduction the model of Lipkin, Meshkov and Glick is a quasi- 
exactly soluble model developed for treating many particle systems. 

The general method for constructing soluble models is based on the incorporation of some 
symmetries of the system which give additional integrals of motion and therefore reduce the 
size of the Hamiltonian matrix to be diagonalized. The Hamiltonian of a many-particle 
system interacting via a two-body force is a sum of linear and quadratic terms in bilinear 
products of creation and annihilation operators related to the quantum states of these par- 
ticles. One starts from the observation that bilinear products of creation and annihilation 
operators can be considered as elements of a Lie algebra, here related in particular to the 
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SU(2) group of the so-called quasi-spin. Lipkin, Meshkov and Glick construct a two iV-fold 
degenerate level model where N is the number of fermions in the system. The two levels are 
separated by an energy e. Here we discuss the simplified version of the LMG model where 
the interaction contains only terms which mixes configurations. In the following we use the 
notations of Ref. 0. Accordingly we introduce fermion operators /3+,/3 m that create and 
annihilate holes in the lower level and a^, a m (m = 1,2, ...,N) that create and annihilate 
particles in the upper level. These operators satisfy the anticommutation relations 

{«m , a+,} = {(3 m , /?+,} = 5 mm , , (1) 

as well as the commutation relations 

k > Pm>] = [dm , = i a m > Pm>] = [<*+ , = . (2) 

The ground state |0) is defined by 

oi m |0) = P m |0> = 0. (3) 

Then the bilinear products 

1 1 N 

j = --N+- £ (a+ « m + /?+/3 m ) , (4) 

m=l 
N 

U =Yj a mPm , ( 5 ) 
m=l 

JV 

J- = a m P™ , (6) 

m=l 

form an su(2) algebra. The Hamiltonian under study can be written in terms of these 
generators as 

H LMG = e 3o + V UI+3-) , (7) 

where e is the separation energy between the two levels, as introduced above, and V is the 
interaction strength. For later purposes it is convenient to introduce the strength parameter 
S instead of V. They are related by 



V 



S e 



(8) 
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The invariant operator of the su(2) algebra 



j 2 = l/2(j + j_+j_j + )+j 2 



(9) 



commutes with the Hamiltonian and provides a constant of motion. Thus the Hamiltonian 
matrix breaks up into submatrices each associated with a different value of j and of order 
2j + 1. Each state in a j multiplet has a different number of excited particle-hole pairs. The 



words it can only change the eigenvalue of jo by two units. The eigenstates of Hlmg have 
therefore an important property, namely their structure is compatible with Hartree Fock 
solutions, so one can study the validity of the Random Phase Approximation, often used 
in the treatment of a system of fermions, against the exact solution. From the definitions 
(HD~(i) it follows that the eigenvalues of jo are given by half the difference between the 
number of particles in the upper level and the number of particles in the lower level. Then 
the maximum eigenvalue of jo an d of j is N/2. Thus LMG conclude that the largest matrix 
to be diagonalized is of dimension 2j + 1 = N + 1. 

Then fixing the number of particles N, LMG diagonalize the largest Hamiltonian matrix 
associated with j = N/2 for several cases. For N = 2, 3, 4, 6 and 8 analytical solutions are 
provided. In addition, the eigenvalues of the multiplet j = N/2 are found numerically for 
iV = 14, 30 and 50. Here by using the polynomial algebra technique we extend the study of 
the LMG model to its entire spectrum. We show that for a given number of particles there 
are two types of states: 1) states with j = N/2, for which the interaction entirely lifts the 
degeneracy. One of these states corresponds to the lowest eigenvalue. These are the states 
analysed by LMG and they belong to the largest matrix to be diagonalized of dimension N 
+ 1 ; 2) states with j < N/2, for which the eigenvalues are identical with those of a system 
with N - 2, N - 4,... The only difference is that these states are degenerate in a system with 



interaction (|7|) mixes the states within the same j multiplet but cannot mix states having 
different eigenvalues of j 2 . It can only excite or de-excite two particle-hole pairs or in other 
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N particles but not degenerate in a system with N -2, N - 4, ...particles, because there they 
are states of type 1. Thus finding the eigenstates of a system with iV particles reduces to 
the diagonalization of the largest matrix once the states of the system with N - 2 particles 
is known. 

In the context of the deformed polynomial algebra we show that the largest matrix 
associated to a given iV can be split into two submatrices of dimensions N/2 + 1 and N/2 
for iV even and two submatrices, both of dimensions (JV+l)/2 for iV odd. These submatrices 
correspond to specific values of the Casimir operator of the deformed algebra. The same 
statement holds for the largest matrix of a system of N - 2 particles, and so on. Alternatively, 
for any j multiplet the corresponding Hamiltonian matrix can be split into two submatrices 
irrespective of the number of particles. These findings are illustrated in detail for the cases 
of iV = 7 and N = 8 particles. 

The splitting of a j multiplet into two submatrices is entirely consistent with the property 
of the LMG interaction (0) that it can excite or de-excite only two particle-hole pairs. 
Accordingly, for A" even, on the one hand the states with 0, 2, N — 2, N particle-hole form 
an N/2 + 1 dimension invariant subspace and on the other hand the states with 1, 3, AT — 1 
particle- hole excitations form another N/2 dimensional invariant subspace. For A" odd the 
states with 0, 2, N — 1 particle-hole excitations and the states with 1, 3, A" particle-hole 
excitations form two distinct invariant subspaces both of dimension (N+l)/2. The deformed 
polynomial algebra provides a "quantum number" denoted here by c to distinguish between 
the eigenvalue of fl7|) for A" even and A" odd. For A" even one has c = and for A" odd 
c = ±1/4. 

Moreover the polynomial algebra technique leads to new representations corresponding 
to new eigenvalues. Some of these are appropriate to a generalized type of LMG model, 
some others are meaningless (see Sec. V). 

As far as physics is concerned one should mention that the LMG model posesses states 
of collective excitations related for example in nuclear physics to giant resonances. For this 
reason Lipkin Meshkov and Glick studied cases with a number of particles and interaction 
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strength relevant to the treatment of nuclei by the Random Phase Approximation. In further 
studies they tested the method of linearizing the equations of motion and the diagram 
summation approximations against exact solutions of their model given in ||. Recently a 
more general Hamiltonian which can test two types of elementary excitations instead of one 
has been proposed by Lipkin ||. 

III. DEFORMED POLYNOMIAL ALGEBRA 

Here we present the LMG model in the context of the deformed algebra. In doing so we 
follow the general method based on the polynomial deformations of the Lie algebra sl(2, R) 
as developed in Ref. 0. 

We start by noting that the Hamiltonian (|7|) is a particular realization of the more general 
Hamiltonian 

H = e [2J + 5(J + + J_)] , (10) 
with the parameter e as defined above, 5 defined by Eq. (H) and 

■Jo = \ Jo, J± = 2^ 3± ■ (11) 
One can show that the operators (|TTD satisfy the following algebra 

[Jo,J±] = ±J±, (12) 

16 2 

[J + , J_] =-—J* + —(2f + 2j - 1) Jo , (13) 

where j is an eigenvalue of @. These relations define a particular case of a deformed 
polynomial algebra as studied in Ref. with the polynom in Jq in the right hand side of 
(Pf) having the coefficients (see Appendix) 

16 2 
« = -^> ^ = °> 7 = ^(2j 2 + 2j-l), A = 0. (14) 

The Casimir operator of this algebra is given by 



C-J J 4 7 4 1 8 7 3 1 2j 2 + 2j-5 2j 2 + 2j-1 

The algebra ([12])- (0) has two types of representations relevant for our discussion. They 
are labelled by q = 1 and q = 2 respectively. More precisely, the q = 1 representations are 
defined by the equations 

J \ JM) = (M + c) \ JM) , 
J+ \ JM)=f(M) \ J,M + 1) , 

J- \JM) = g(M) \J,M -1} , (16) 
with M = — J, J, J = 0, -, ... , c e 3ft and 

f(M - 1) «/(M) = ij(J - M + 1)(J + M) 

x [2j 2 + 2j - 1 - 4J 2 - 4J - 4M 2 + 4M + 8(1 - 2M) c - 24c 2 ] , (17) 



where the real number c, constrained by Eq. ( fj2|) (see Appendix), can take three distinct 
values given by 

c = 0, (18) 

or 

c = ±[^(j + l)-^-^+l)] 1/2 - (19) 
The q = 2 representations are defined in an invariant subspace satisfying 

J \J'M') = ^- I J'M') , 



J+ I J'M') = /'(M') |J',M' + 2) , 
J_ |J'M') = #'(M') |J',M' - 2) , (20) 



where J =0,1,2,... and 



f'(M' - 2)g'(M') = -L^/ - M' + 2)(J' + M' 



x (2j 2 + 2j - 1 - J 2 - 2J - M + 2M ) , (21) 



if M = - J , J -2, J and 



f'(M' - 2)g'{M') = _L(./ -M' + 1)(J' + M' — 1) 



x (2j 2 + 2j- J' 2 -M' 2 + 2M') , (22) 



if M = -J + 1, J -3, J - 1. 
/ 1 3 

The cases J = are particular in the sense that J' must be equal to j. The 

relations satisfied by the basis vectors \jm) are 

\jm) = —\jm) , 
J+ = /O) |j,m + 2) , 

J- \jm) = g'{m) \j, m - 2} , (23) 

with 

/'(m - 2)^'(m) = -^(j + m)(j + m - l)(j - m + l)(j - m + 2) . (24) 

The Hamiltonian ([?[) can be associated to the representation q = 2 since in this case the 
invariant subspace is spanned by the vectors \jm) on which the deformed generators act as 
follows 

MJ m ) = , ( 25 ) 

J + \jm) = -^—sj (j - m - l)(j - m)(j + m + l)(j + m + 2) |j,m + 2) , (26) 



J_|j m ) = _ + m - l)(j + m)(j - m + l)(j - m + 2) |j,m - 2) . (27) 

One can see that these relations can be recovered from the equations ( p3|) if j is a half integer 
and from fl2"0|) but with J = j if j is an integer. 

If one now calculates the eigenvalues of the Casimir operator ( |1"5|) for the representation 
q = 2 and J = n, i.e. an integer, one gets 



S 



(C) 



j' =n,q=2 



1 ^^ + 2)[j(j + l)-^ + l) 2 ] ; 
(" "IK"" 2)[i0' + l)-^ 2 ] ; 



M = n,n — 2 — n 



M = n — 1, n — 3, — n + 1 



(28) 



One can therefore see that C has two distinct eigenvalues in the space spanned by | J M >. 
This shows that the representation q = 2 is reducible. One can easily prove that it can be 
split into the direct sum 



(/ = n,q = 2) c=0 



2 c =o 2 c =o 



(29) 



i.e. the q = 2 representation can be decomposed into two q = 1 representations and this 
takes place for c = only. 

A similar decomposition also holds for half integer j. In this case one has 



0' = n+~,g = 2) 

2 c =o 



2 c=l/4 2 c=-l/4 



(30) 



and the eigenvalue of the Casimir operator is the same for q — 1 and q = 2 

' ./(./ • l.H./ • 2) . 



(31) 



( C )j=n+l/2 4i y 2 

We can then conclude that the Hamiltonian matrix of each j multiplet can be split into two 
submatrices. Examples are shown in the next section. 

^From now on we are concerned with q — 1 representations only. Then in the invariant 
subspace defined by Eqs. (|Tjl) we obtain the following Hamiltonian matrix to be diagonalized 



(H) 



2J + 2c 5f(J-l) 

6g(J) 2J-2 + 2c 5f(J — 2) 

Sg(J-l) 2J-4 + 2c 6f(J-3) 

Sg(J-2) 2J-6 + 2c 










-2J + 2 + 2c Sf{-J) 
5g{-J+l) -2J + 2c 
(32) 
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The diagonalization amounts to solving the secular equation 

det\(H) ij -E5 tJ \ = (33) 

IV. THE CASES N = 7 AND 8 PARTICLES 

In this section we give analytic and numerical results for the eigenvalues of (0) obtained 
by solving the eigenvalue equation (|33"D for 7 and 8 particles. For a larger number of particles 
we checked that we perfectly agree with the numerical values of Ref. 0. 

A. N = 7 

For N = 7 there are 2 7 =128 states. The largest Hamiltonian matrix corresponds to 
j = 7/2. According to the relation ( |30| ) where we have to take n = (N — l)/2 this matrix 
splits into two equal matrices, both having J = (N — l)/4. The same procedure applies to 
the j — 1 multiplet which is the largest multiplet of N = 5 particles, and so on. In Table 1 
we give the possible j multiplets, their multiplicities rrij and the corresponding values of J. 
Analytic forms of the eigenvalues can be easily obtained only for j = 1/2 (trivial case) and 
j = 3/2. For j = 5/2 they are numerically obtained from the secular equation (|33D which 
in this case becomes 

- 6cE> - (H + *#)E + H a + = (34) 

with c = ±1/4. For j = 7/2 the secular equation (|33D leads to 

„a /17 18 ^o. ,_,r, . 888 ^r,, „ 105 75 135 ~a „ /__\ 

E A - 8cE 3 - — + —5 2 )E 2 + 38 + 5 2 )cE + + — 5 2 + 5 4 = 35 

v 2 7 ; V 49 ; 16 14 343 V ; 

with c = ±1/4 which is also solved numerically. The dependence of the eigenvalues on the 
parameter 5 is exhibited in Fig. 1. 
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B. N 



For N = 8 there are 256 states. This case is more fortunate because one can push the 
analytic calculation further than for N = 7. The largest Hamiltonian matrix corresponds 
to j = N/2 = 4. According to where we have to take n = N/2 this matrix is split into 
two submatrices one corresponding to J = N/4 = 2 and the other to J = (N — 2)/4 = 3/2. 
The next to the largest multiplet, with j = N/2 — 1, can be seen as the largest multiplet 
corresponding to iV — 2 particles and one can use the decomposition (|2D| ) again. In the 
N — 2 particles case all the eigenvalues are nondegenerate while in the TV particle case the 
same result is valid for j = N/2 — 1 but with some degeneracy for the eigenvalues. In Table 
2 we exhibit the multiplets j and their multiplicities rrij, the values of J consistent with 



(P9|) for each j and the corresponding analytical solutions for the eigenvalues obtained from 
the secular equation (|33|) . One can check the consistency of these analytic expressions with 
those of Ref . by using the relation @ which gives 5 in terms of V/e [|J . 

In Fig. 2 we plot all positive eigenvalues of (0) as a function of the strength 5. One 
can notice an expected degeneracy at 5 = and some degeneracy at large values of 5. In 
particular the largest eigenvalue with j = 4, J = 2 becomes degenerate with the largest 
eigenvalue with j = 3, J = 3/2. The same eigenvalues with opposite sign play the role of 
the ground state of the system, which thus becomes degenerate for large value of 5. 



V. SUPPLEMENTARY EIGENVALUES 

We have used the sl(2,R) deformed algebra to study the spectrum of the Hamiltonian 
(0). By construction, this algebra is richer than the su(2) algebra. Its representations have 
three labels instead of one, as for su(2). Thus the number of representations is larger. 
Moreover we can see that once j and c are fixed in (|I7|) , one can choose value of J different 
from j such as the right hand side remains always positive. These values have nothing to 
do with the LMG model. However it turns out that some eigenvalues of Eq. (^) associated 
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to these representations are quite similar to those of ([?]). It would be interesting to find out 
if they can be related to a more general Hamiltonian than LMG. This is the subject of a 
further study. 

VI. SUMMARY 

Here we have presented a calculation of the whole spectrum of the Lipkin-Meshkov-Glick 
Hamiltonian presented in the context of deformed polynomial algebra. For any given number 
of particles iV the spectrum first divides into j multiplets of the su(2) algebra. The eigenval- 
ues associated with the largest j are nondegenerate except for E = 0. We have shown that 
the Hamiltonian matrix of each j further splits into two submatrices corresponding to two 
distinct irreducible representations of the polynomial deformed algebra. These representa- 
tions bring new "quantum numbers" , one of them allowing us to distinguish between N even 
and N odd. In order to illustrate the method we have derived explicit analytic expressions 
for the eigenvalues of the LMG Hamiltonian for N = 7 and 8. Our method can be extended 
to any N. 

Furthermore, we have shown that the deformed polynomial algebra related to the LMG 
model implies a larger spectrum than that of the model itself. Some of the new eigenvalues 
present characteristics similar to those of the LMG model itself and may lead to a kind of 
generalized model. 

We hope this study could shed a new light into the LMG model and could inspire further 
applications. 

Appendix 

For the purpose of selfconsistency, here we recall a few results obtained in Ref. which 
are directly exploited in Sec. 3. One of the aims of was to construct finite dimensional 
representations of the polynomial deformed algebra 

[Jo, J±] = ± J± (36) 
12 



[J+, J_] = a J 3 + J 2 + 7^0 + A , a,/3, 7 ,A e 3? (37) 

related to quasi-exactly soluble models. Such representations imply the existence of kets, 
shortly denoted by \JM), such that 

Jo \JM) = (M/q + c) \JM) , 
J + \JM) = f(M) \ J,M + q) , 
J_ I JM) = g(M) \J,M-q) (38) 

In fact the kets are characterized by four labels: the number J = ^ 1, ... related to an 
eigenvalue of the Casimir operator 

C = J + J. + Hj} + (|-H)JS + (£-| + I)JJ + (|-l + A)Jl, (39) 

which gives the dimension 2 J + 1 of the representation, the number M = -J, —J + 1, J 
which represents the eigenvalues of Jo, the positive integer q connected to the strength of 
the raising and lowering operators J± and the shift real number c which also enters the 
eigenvalues of C . Thus c depends on J. 

The highest weight vectors of each representation impose the constraints 

f(J) = f(J-l) = ... = f(J-q + l) = 0, 
g(-J) = g(-J + 1) = ... = g(—J + q - 1) = . (40) 

These relations ensure that the dimension of the representation is 2 J + 1 and they lead to 
the following system of q equations for the number c 



.3 



3(d-l) 2 J 2 - J(d + l) + l 2 -dl + d 2 2J-d-l 

— -c z H — c H < 

2q q 2 2q 



2J 2 (d -I)- 2J(d 2 - I 2 ) + d 3 - d 2 l + dl 2 - I 3 l 2 -d 2 + 2J(d-l) ] 

+ 4q~ 3 + 4q~ 2 1 

n[2 (d-l)c J 2 - J(d+l)+d 2 -dl + l 2 2J-d-l), 
+P [c 2 + + ^ !—z + 



3q 2 6q 

+7(c+^)+A = (41) 
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with I = 0, 1, q — 1. The non-negative integer d, introduced above for convenience, has to 
take specific values according to J and I. These values have been given in J7[. For the real 
value of a, (3, 7 and A of Eqs. ([TJ]) the obtained system of equations is compatible if and 
only if q = 1 or q = 2. 

For q = 1 one has d = and one remains with a single equation 

ac [c 2 + J(J + 1)] + /3 [c 2 + J(J+1) ] + 7 c + A = 0. (42) 

For q = 2, one has: 1) d = when / = and J is an integer or when / = 1 and J is a half 
integer, 2) d — 1 when I = and J is a half integer or when I = 1 and J is an integer, with 
corresponding equations of type ([41]). 

The representations (|38|) are completely specified with 

f(J-kq-q-l)g(J-kq-l) = (k + l) 

{a (til + C )3 + ( J^ + c) 2 +7 + c) + A 

g g g 

-i[3a (— - + c) 2 + 2(3 (— - + c) + 7] fc 

2 g g 

+ ^[3a (— ^ + c)+/3] fc(2A; + l) 
g 

-^ 2 (fc + l)} (43) 
where I = 0, 1, g — 1 and fc = 0, 1, (2 J — d — l)/q. 
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TABLE I. Eigenvalues of the Hamiltonian (0) for N = 7 particles. 
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TABLE II. Eigenvalues of the Hamiltonian (^) for N = 8 particles. 
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FIG. 1. Positive eigenvalues of the LMG Hamiltonian ([?]), as a function of the parameter 5 
defined by Eq. (||) for N = 7. The eigenvalues are labelled by j, J and correspond to the rows 3 
and 4 of Table 1. 
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FIG. 2. Positive eigenvalues of the LMG Hamiltonian ([<]), as a function of the parameter 5 
defined by Eq. (||) for N = 8. The eigenvalues are labelled either by j, J or by j, J, sign when 
necessary, where sign means the sign in front of the inner square root in the last column of Table 
2. 



18 



